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In this study we characterize the rheology of fluidized granular matter subject to secondary forcing. 
Our approach consists of first fluidizing granular matter in a drum half filled with grains via simple 
rotation, and then superimposing oscillatory shear perpendicular to the downhill flow direction. 
The response of the system is mostly linear, with a phase lag between the grain motion and the 
oscillatory forcing. The rheology of the system can be well characterize by the GDR-Midi model 
if the system is forced with slow oscillations. The model breaks down when the forcing timescale 
becomes comparable to characteristic time for energy dissipation in the flow. 
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Describing and accurately predicting the flow of granu- 
lar matter is important for many industrial applications, 
as well as geological processes [IHl]- For example, ac- 
curate models for the flow of dense granular matter are 
important for applications such as soil compaction [S], 
rock avalanches [SI [7], or manufacturing of pills from 
powders fS^. However, unlike fluids where flow is well 
described by the Navier-Stokcs equations, no such set of 
equations has been shown to predict the large variety 
of observed granular flows. Nevertheless granular flows 
tend to be very reproducible and describeable by simple 
empirical continuum equations. 

Recently significant progress has been made to classify 
a wide range of granular flow geometries within a unify- 
ing empirical framework, da Cruz et al. and the GDR- 
Midi group have developed a model employing a single 
dimensionless parameter, the inertia number, that char- 
acterizes the flow behavior based on the importance of 
inertia relative to the confining pressures [SHU] ■ Key as- 
pects of the GDR-MiDi model are the presence of a yield 
stress and the approach of an effective viscosity from a 
non-zero constant towards a larger asymptotic value in 
the limit of large strain rate. The GDR-Midi model has 
been successfully used to develop constitutive equations 
that accurately predict shear fiow profiles for various ge- 
ometries such as chute flow and rotating drum flow. 

Test of this model and other models of granular 
flow [13] have focused on 2D flow profiles. In such flows, 
from a microscopic view of forces transmitted through 
particle contacts, the fabric tensor and force chains [TJ] 
maintain direction or rotate within a plane. Since the 
GDR-Midi model was developed with such flow profiles it 
is likely that this effect is contained in the expression for 
the effective viscosity. However, many relevant bound- 
ary conditions in granular flows involve additional forces 
perpendicular to the 2D flow proflle and 3D flow fields. 
Common examples include rock avalanches flowing down 
real surface topographies, or grains exiting the orifice of a 
hopper. From the microscopic view the fabric tensor and 
force chains are now forced in 3D rather than a plane. 

In this paper, we present a new table top geometry 
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FIG. 1. (Color online.) Schematic of the experimental appa- 
ratus. The drum is half filled with glass beads of diameter 
d and density p, and the grains are fluidized by rotating the 
drum with angular frequency Q. A secondary shear is applied 
by oscillating the drum with frequency /, and the oscillating 
motion of the grains at the surface is recorded with a high 
speed camera. The plot to the right shows the motion of the 
drum (dark, blue) and 2mm grains (light, red) measured for 
/ = 2.5 Hz; plus sinusoidal fits (solid lines). 



to characterize the rheology of granular matter that is 
already exhibiting a shear flow, with the aim to eluci- 
date how flowing granular matter responds to additional 
forcing. Our approach consists of first fluidizing granu- 
lar matter by rotation in a horizontal cylinder rotated 
about the cylindrical axis, as shown in the schematic in 
Fig. [l] In addition to rotation, shear is applied orthogo- 
nal to the flow direction via horizontal oscillation of the 
entire cylinder. To test the GDR-Midi model we com- 
pare our results to the prediction of the model. We find 
that the model is unable to predict the rheology at large, 
fast forcing, but is able to reproduce the correct scaling 
relationship and predict the rhoelogy at low forcing. 

In our experiment, we rotate a drum of diameter 10 
cm and length 67 cm half full of monodisperse beads 
about the long axis at frequency fl. We use glass beads 
(Glenn Mills, Inc.) with diameter d. A high speed cam- 
era (Photron FastCam) placed above the drum records 
the motion of the grains at a frame rate of 500 frames/sec 
for 2.096 seconds. The camera is located 10 cm from the 
end of the cylinder to avoid boundary effects from the 
end caps. 
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FIG. 2. (Color online.) (a) Phase lag 5(j!)(/) vs oscillation 
frequency for representative particle sizes d and rotation rates 
fj. The units for d and fi are mm and Hz, respectively. 



To assess the rheology of the flowing layer a secondary 
forcing is superimposed by driving the drum sinusoidally 
along the long axis at driving frequency / (see Fig. IT]). 
We define the x-axis to be parallel to the oscillating mo- 
tion and the z-axis to be perpendicular to the granular 
surface. 

Using images taken with the high speed camera both 
the motion of the grains and the drum are tracked si- 
multaneously using an algorithm by Crocker et al. |15j . 
The black stripes seen in the middle of the image is tape 
placed at the surface of the drum to identify the drum's 
horizontal position. Figure [T] shows a plot of both the 
horizontal x-axis motion of grains at the surface and the 
horizontal motion of the drum, respectively. 

The velocity of the drum and the grains along the drum 
axis in the x-direction u^ (see Fig. II]) are fit to a sinu- 
soidal function w^. = —2tt A f sm{2TT ft + 0) to obtain a 
frequency /, amplitude A, and phase of the particle 
and drum motion. The phase of particle motion (j)p and 
phase of drum motion (f>]j differ. The peak acceleration 
due to horizontal shaking is less than 0.02g, significantly 
smaller than the downhill acceleration of particles and 
too small to fluidize the material directly when the tum- 
bler is not rotated. However, for our larger oscillation fre- 
quencies near 4 Hz the downhill velocity becomes slightly 
perturbed. 

The frequency (/) dependent response of the fluidized 
granular matter to oscillatory forcing can be described 
with the phase lag Scj) — (j)£, — (pp. Figure l2ja) shows the 
frequency dependence for a few representative values of d 
and fl. The data shows that the phase lag varies linearly 
with oscillation frequency / and goes to zero as / — >■ 
within experimental error. 

We note that an oscillating layer of Newtonian fluid 
would exhibit a scaling of Scj) ^ y/f, unlike the linear 
scaling we observe in our data |16j . This discrepancy 
implies that in order to accurately model the behavior of 



this system a non-newtonian model is needed. 

Thus we apply the GDR-Midi model that has been suc- 
cessfully used to model simple tumbler flows. The GDR- 
Midi model is a continuum model treating the granular 
medium as a frictional material with a shear rate and 
pressure dependent friction coefficient. We simulate the 
model as layers of solid planes slipping past each other 
and interacting through friction, similar to a da Vinci 
fluid model [T^ . When there is a velocity difference be- 
tween adjacent layers they exhibit stresses on each other 
a = /iPsign(7), where /x is a strain rate dependent fric- 
tion coefhcient, P is the pressure, and 7 is the strain rate. 
The strain rate dependence of fi is written as 



fj- = Us + ifJ-2 - M.s)/(1 + Iq/I), 



(1) 



where /ig, ^2, and Iq are material dependent constants, 
and the inertia number / = \j\d/ y/P/p, where d is the 
particle diameter and p is the particle density. This 
model can be recast into a Non-Newtonian fluid model 
by defining an effective viscosity as rj — p{I)P/\'^\ such 
that a = Tj'^. 

To reduce the complexity of our system we will treat 
the rotating drum apparatus as a oscillating inclined 
plane. In this geometry, the GDR-Midi model can be 
written as the following set of coupled differential equa- 
tions 
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where the notation Uij is the partial derivation of the ve- 
locity in the i direction with respects to the j direction, 
7 is the strain tensor, and 9 is the angle of inclination of 
the surface. The first equation governs the flow down the 
inclined plane and the second equation governs the flow 
in the oscillatory direction. In the absence of secondary 
forcing (typical incline plane flow) Jop et al. [TH] showed 
that p{z) is constant with depth in the steady state. As a 
consequence, the inertia number is also a constant inde- 
pendent of z where / = [(/i2— /Xs)/(/x— /is) — l]~^/o. Using 
the definition of inertia number and P — pgzcos{9), one 
can invert the previous equation to find 7(2) — C\fz^ 
where the constant C = I/{dy/gcos{9)). Therefore, in 
the absence of secondary forcing there exist an analytical 
form describing the flow. For the fully coupled equations 
shown in Eqn. [2] there is no analytical solution. Instead 
we must use numeric integration to solve for the flow pro- 
file, and we have checked that our numerical integration 
reduces to the correct solution in the absence of oscilla- 
tory flow. 

Before we continue to the results of the model we seek 
to understand if there is a natural time scale t^ present 



in the system so that we can compare our forcing time 
scale tf = 1// to the relaxation time of the material. To 
find t^ we consider the differential equation describing 
the downhill flow in the absence of secondary forcing. In 
this differential equation 77 can be written explicitly as a 
function of z using our results from before, where rj — 
fipgcos{9)y/z/C. Making a change of variable z' = z/h, 
where h is the thickness of the flowing layer, and plugging 
T] into the differential equation for the downhill flow we 
find 



duy{z') fj.gcos{9) d 



dt 



Ch^/^ dz- 



- [yz'uy^^.'j +gsin{e). (3) 



By dimensional arguments the above equation has a char- 
acteristic time scale t^ = Ch^^^/{fj,gcos{6)), where i^ is 
a characteristic time scale for the flowing layer to dis- 
sipate it's energy. For instance, if we consider the case 
of a flowing granular layer of thickness h at an inclina- 
tion angle 9, i,, is a characteristic time for the grains 
to stop flowing if we suddenly changed the inclination 
angle to 9 — 0°. We also note that by integrating the 
strain rate, the downhill velocity at the top of the flow- 



ing layer is Uy{z = 0) 



Uto 



2/3C/i3/2 and therefore 



Ui = 3utop/{2figcos{9)). 

Using the two time scales t^ and tf we can define a 
dimensionless number tri/tf which characterizes the in- 
fluence of the secondary forcing. When t^/tf <C 1 the 
secondary forcing takes place over a much longer time 
scale than the relaxation time of the grains and therefore 
during each oscillation the grains can relax. However, 
when t,f/tf > 1 the secondary forcing is fast enough that 
during each oscillation the grains will not be able to com- 
pletely relax. We may expect that as t^/tf approaches 
unity the response of the material may change. 

Now knowing the dimensionless parameter describing 
the system we numerically integrate the coupled differen- 
tial equations Eqn.[2]to solve for over 800 different S(j){f) 
response curves for various combinations of parameters 
spanning the ranges: d = 1 — 3 mm, 9 = 18° — 28°, 
h^7-15 mm, /i^ = 0.33 - 0.42, ^2 = 0.55 - 0.64 and 
lo — 0.279. The results of the model are shown in Fig. [3] 
as the black solid line, where the response curves are 
plotted as a function of trj/tf. The black line shows that 
all 800 response curves collapse to a single curve when 
rescaled by the dimensionless ratio tjj/tf, indicating that 
tri/tf is an appropriate measure to characterize the in- 
fluence of the secondary forcing. The data have been 
plotted on a log-log scale to show that at small tr,/tf the 
model predicts a linear response regime. 

To compare our experimental results to the model we 
must compute t,i/tf — 3utopf / i^ng cos{9)) for our exper- 
imental data. However, for our experiments we do not 
know the value of /i, though we do know the other values 
^7 5: Utop: smd f. Therefore we treat /x as a fitting param- 
eter and find /i for each experiment by fitting 5(j}{f) to 
the linear response regime found from the model. The ex- 





,' 

• 
• 


10° 


^^^ 




i^ 


ff- 


* d,Q= 1,0.09 , 




>^4" 


■ d,Q= 1,0.13 




A* I 




• d,Q=l,0.18 


in-' 


- xl 






♦ d,Q= 1,0.28 , 




X ' 






• d,Q = 2,0.18 : 




y . . 




1 


• d,Q = 3,0.18 : 



10" 



10"' 

trtjtj 



10" 



FIG. 3. (Color online.) Dimensionless response curves are 
shown for our model and the experimental data. The black 
solid line is data from the model and the black dash line is 
an extension of the linear regime found in the model and is 
meant as a guide to the eye. The symbols are experimental 
data. The units for d and fi are mm and Hz, respectively. 



perimental results are shown as symbols in Fig. [31 where 
each data set was fitted to the black dash line to find 
a value for /i. The fitted values of [i ranged from 0.45 
- 0.54. In the GDR-Midi model for steady state incline 
flow [I = tan(0) [18]. For our experiments, 9 ^ 23° giving 
an expected value of /i ^ 0.43 close to the values found 
for our experiments. 

Figure [3] shows that the linear response regime of the 
experiment extends beyond the linear response regime 
predicted by the GDR-Midi model. Since the experi- 
mental data at larger trijtj have a larger phase lag than 
predicted, it suggest that the effective viscosity of a real 
granular material upon large secondary forcing is smaller 
than the model predicts. One possibile reason why the 
material exhibits less flow resistance under fast forcing 
may be that granular flows are transliently weak when 
their principal stress axis is changed [19!, likely since fi- 
nite strain is needed to realign the fabric tensor and force 
chains in a new direction. Our experimental results in- 
dicate that the material remains weak up to a timescale 
comparable to the visous relaxation timescale. 

With our new method now established, we also inves- 
tigate the behavior of two binary mixtures of 1 mm and 
2 mm diameter particles. Figure [4] shows the response 
curve for two different mixture ratios of 1 mm to 2 mm 
beads at 17 = 0.18 Hz. To ensure that the grains had 
fully segregated we allowed the drum to rotate at least 
50 times before recording the flow. 

The curves in Fig. |4] show that the mixtures behave in 
a manner that is bound between the behavior observed 
for pure 1 mm or 2 mm bead systems, and exhibit a non- 
linear response with increasing / which is much different 
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FIG. 4. (Color online.) (50 is shown for two different binary 
mixtures of 1 mm and 2 mm particles. The blue (dark gray) 
data points are for a mixture of 75% 1 mm beads by volume 
and the orange (light gray) data points are for a mixture of 
50% 1 mm by volume. Data for a the monodisperse systems 
are also shown. 



than the monodisperse system and the GDR-Midi model. 
We find that for small / both systems behave like a sys- 
tem of pure 1 mm beads. This is surprising given that 
after segregation very few 1 mm particles will be in the 
fluidized layer, but yet these few 1 mm particles can sig- 
nificantly weaken a 2 mm sample. However, for higher 
/ we see that the system containing 50% 1 mm beads 
converges to the pure 2 mm bead curve, while the sys- 
tem containing 75% 1 mm particles has a response that 
is between the behavior observed for the 1 mm and 2 
mm monodisperse systems. This suggests that with in- 
creasing oscilation frequency more 1 mm beads may be 
involved in the oscillatory flow field. 

In this paper we have introduced a new approach to 
quantify the response of fluidized grains to secondary 
forcing. Our new apparatus consist of first fiuidiz- 
ing granular matter by rotation in a horizontal cylin- 
der rotated about the cylindrical axis, as shown in the 
schematic in Fig. [l] In addition to rotation, shear is ap- 
plied orthogonal to the fluidizing mechanism via horizon- 
tal oscillation of the entire cylinder to probe the rheologi- 
cal response of the fluidized grains. We find that fluidized 
monodisperse spheres at low forcing are well described 
by the GDR-Midi model and that the response function 
S(j>{f) can be used to find the effective friction coefficient 
fj, for steady state shear. However, at faster forcing fre- 
quencies the GDR-Midi model breaks down, predicting 
a larger effective viscosity than observed in the experi- 



ment. We attribute the smaller than expected viscosity 
observed in experiments to the slow adaptation of the 
fabric tensor and force chain directions to the direction 
of forcing. Unlike traditional rheometers that elucidate 
materials properties from the amplitude and phase of the 
response to oscillatory forcing, we suggest that the rhe- 
ology of the flowing state is best determined via small 
perturbations of a primary flow field. Our device thus 
provides a template for a simple rheometer for the fiow- 
ing state of granular matter. 
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